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An operator form of the 3N bound state is proposed. It consists of eight operators formed out of 
scalar products in relative momentum and spin vectors, which are applied on a pure 3N spin 1/2 
state. Each of the operators is associated with a scalar function depending only on the magnitudes 
of the two relative momenta and the angle between them. The connection between the standard 
partial wave decomposition of the 3N bound state and the operator form is established, and the 
decomposition of these scalar function in terms of partial wave components and analytically known 
auxiliary functions is given. That newly established operator form of the 3N bound state exhibits 
the dominant angular and spin dependence analytically. The scalar functions are tabulated and 
can be downloaded. As an application the spin dependent nucleon momentum distribution in a 
polarized 3N bound state is calculated to illustrate the use of the new form of the 3N bound state. 
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I. INTRODUCTION 



Today it can be considered as standard to solve the Schrodinger equation for three nucleons numerically with high 
precision. This can be done either in the form of Faddeev equations |lj,|2|, usm S hyperspherical expansion |3j.|4( or the 
Gaussian-basis method 5]. Also the Greens- Function-Monte-Carlo (GFMC) and No-Core-Shell Model approaches 
have been applied to the three nucleon system [|| Q. These methods except for GFMC have always been based 
on standard partial wave expansions. Only recently for the study of three bound bosons the Faddeev equation has 
been solved directly in terms of relative momentum vectors 0, 0, E3 ■ Though any expectation value for a three- 
nucleon (3N) bound state can be determined numerically from wave functions obtained by any of the above mentioned 
methods, no analytic insight can be extracted for the spin and momentum (or configuration) dependence of the 3N 
bound state. 

For the deuteron Rarita and Schwinger [Til IT^ introduced an operator form of the deuteron state in terms of spin 
operators and the relative position vector. This representation is ideal to exhibit the probabilities for finding any spin 
orientations in relation with the relative position (or momentum) vectors of the two nucleon in a polarized deuteron. 
In the case of a momentum representation of the deuteron derived from modern nucleon-nucleon (NN) forces this is 
illustrated in Ref. |r| and for a coordinate representation in Ref. [Tij |. 

For three nucleons it is much more difficult to express the spin and momentum dependence of the bound state in an 
analytic form. Again, this analytic form has been worked out long time ago by Gerjuoy and Schwinger |l5j |. For the 
3N bound state, however, the functions multiplying the different scalars built from spin and position vectors depend 
on three variables, namely the magnitudes of the Jacobi vectors and the angle between them. In Ref. 01 t ne wave 
function is expanded in nine operators with nine corresponding scalar functions, which were numerically unknown in 
those days. It is the aim of this paper to establish an analytic link between these functions and the usual partial wave 
representation of the 3N bound state. This will lead to an alternative representation of the 3 H ( 3 He) bound states, 
which is more accessible to analytic insights into the spin structure of the 3N bound state. In our analysis we also 
find that the last term given in Ref. jl5| is redundant, further simplifying this representation. 

The paper is organized as follows. In Section II we start from the standard partial wave representation of the 3N 
bound state and reformulate it in terms of scalar operators acting on pure spin states with the correct spin quantum 
numbers of 3 H ( 3 He). This reformulation then leads at the same time to the scalar functions. Section III is devoted 
to the numerical investigation of the scalar functions. In Section IV we give an example for the application of the 
operator form of the 3 He wave function and evaluate the probability to find a neutron with given momentum vector 
and polarized in the direction of the overall polarization of 3 He. Several appendices include technical steps for the 
derivations. Finally we summarize in Section V. 
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II. THE OPERATOR FORM OF THE 3 H ( 3 HE) BOUND STATE 

A. Derivation 



The starting point for the derivation of an operator form is the standard partial wave representation of a 3N bound 
state. Here we do not use the isospin formalism and choose particles 1 and 2 to be neutrons(protons) and particle 
3 to be the proton(neutron). For the derivation we will assume a 3 H bound state, but the generalization to 3 He is 
obvious. In momentum space the 3N bound state with total angular momentum 1/2 and magnetic quantum number 
m can be written as 



* m (p,q) 



(pq|* m > 

E l(4)^s)((4)^5|(pq|* r ' 



sSMs 



E I (shsMs) E C{Ls\,m - M s M s )y^ Ms (pq) *i XLsS (pq) 



sSM s 



EE 

sS IXL 



IXL 



yfximi^s) 



^ixLss{pq), 



(1) 



where p and q are the standard Jacobi momenta , and the p and q stand for the corresponding unit vectors. The 
3N spin state \(s^)SMs) is constructed by coupling the spin state of the two neutrons with total spin s and the spin 
of the proton (i) to the total spin S and magnetic quantum number Ms of the 3N system. Furthermore, A, and 
L are the relative orbital angular momenta of the two neutrons (related to p) , the orbital angular momentum of the 
proton (related to q) and the total orbital angular momentum of the three nucleons. The bracket in the last equation 
is a convenient abbreviation for the LS-coupling employed here. The quantities ^ixLss(pq) represent the partial wave 
components of the 3N bound state. They are for example determined by the solution of the Faddeev equations 0,0- 
Typical numbers for a good representation of \E , ' T1 (p, q) are I < 3, A < 3. The wave function is antisymmetric in 
the two neutrons, which constrains I + s to be even. Note that we solve the Schrodinger equation using the isospin 
formalism. The particle basis wave functions, which enter in Eq. are a combination of total isospin T = 1/2 and 
T = 3/2 wave functions, namely 



IXLsS 



(pq) = Vs ( -)^+3/2 



2 t=l T= 
g ^IXLsS 



Hpq) 



1 t=l T= 
3 ^IXLsS 



'(pq) 



(2) 



where Mt is the third component of the isospin, namely Mt = 1/2 for 3 He and —1/2 for 3 H. Note that the isospin 
of the neutron-neutron or proton-proton two-body subsystem is restricted to t = 1. The overall factor \/3 is a 
consequence of the identity of the nucleons in the isospin formalism and insures the correct normalization. 

The right hand side of Eq. naturally decomposes into four parts according to the total spin S =1/2 or S = 3/2 
and the total orbital angular momentum L = 0, 1, and 2 as 



* m (p,q) = EE 
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+EE 

s 

+E 
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^IXUiiPl) 
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Pi 
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Pi 



1/2 



D 1/2 . 



(3) 



At first we consider the 2 5*i/2 part which can be separated in even and odd terms with respect to / as 



2 Si 
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E 



^l(pq) I (o|)i> 
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The first spin state occurring in Eq. will be denoted as | x m ). It carries the correct spin quantum numbers of the 
3N bound state and is explicitly given as 



^\)\) = ^{xtX2-XiXt) X? 



(•5) 



By introducing the spin operator 



<7(12)=-(<r(l)-<r(2)), (6) 

which is odd under the exchange for particles 1 and 2, one can verify by straightforward calculation that 

<r(12) • <r(3) I (O^m) = -V3 | {l\)\m). (7) 

This leads to the second spin state in Eq. Thus, taking the antisymmetry with respect to nucleons 1 and 2 into 
account, Eq. 10} can be rewritten as 



I 2 S 1/2 ) = Ix^E^^^p q^.^/S 



-L<r(12) • ff (3) | X m ) £ ^^(P ' q) *hoi*(p?)- 



(8) 



Here the standard relation of 3^](pq) to the Legendre polynomial Pi(p ■ q) has been used. These are the first two 
examples for the operator form of the 3 H bound state. The state | x m ) with its correct spin quantum numbers for 
3 H is multiplied by scalar functions and occurs either by itself or is acted upon by a rotational invariant expression 
formed out of spin operators. Below additional rotational invariant expressions will appear in the 3 H wave function, 
which are formed out of spin operators and momentum vectors. 
The second part of Eq. J3J, 2 P\/2, has the following explicit form 



^(pq) I (05)5) 



IX 

I 2 P l/2) ls=0 + I P1/2) l*=l 



IX 



*;aiu(w) 



(9) 



First let us consider | 2 Py 2 ) l«=Oi insert the Clebsch-Gordon coefficients and use the standard descending spin 
operator, which is related to the spherical component cr_(3)0l- Doing this one obtains 



2 Py- h ) | s =o= ^= £ (-a (3)^°(pq)+a_(3)^ 1 (pq)) 



x m -*)^u W dpq)- 



(10) 



In the last equation we also used the fact that 3 H has positive parity, which enforces I = A. Next we need to face the 
problem of the infinite sum over I, which includes the angular dependence on p and q. As shown in the Appendix A, 
the following relation holds 

yr(pq)=Q(p-q)^r(pq), (11) 
where the coefficients c;(p • q) are analytically known functions. Thus, we obtain for Eq. I|l()|l the expression 



2p m 



1/2 



n/3 



ao(3)y i 1 1 °(pq)+^(3)X 1 1 (pq)) \X T - 

£ c ;(p-q)*nioite)- 



(12) 



Further one recognizes that 3^i" l (pq) can be expressed in terms of the spherical components of the cross product 
p x q, which turns Eq. l|12f) into the operator form 



1 2 pm \ I 



1 /31cr(3)-pxq 



pq 



£ ci(f>- q) *hio*(p9)- 



(13) 
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The other P-component of Eq. @, | 2 Py 2 ) U=i is a little more complicated. Inserting the Clebsch-Gordon 
coefficients and making use of the relation in Eq. (0 gives 



1/2 



-cr {12) ■ er(3) 



X 5>Hiij(P?) (a (3)^°(pq) - M3)3#(M)) I X m=i )- 

lodd 



(14) 



Using again the relation given in Eq. and expressing the quantity y\™ (pq) in terms of the spherical components 
of p x q leads to the intermediate result 



2 pm 



1/2/ 



1 1 



-In 



-±= £ Q ( P • q) * Kn i(pg)^(12) • ^(3)^-^5 | \ 



m 



(15) 



and finally to 



2 pm \ 
-^1/2/ 



1 1 

47T 

1 



x-Kl2).pxq-(ff(3)x<r(12))-(pxq)) | x m ). (16) 
The next term to calculate from Eq. © is | PJjo)- Again we insert the Clebsch-Gordon coefficients and find 



4 r,m=l/2\ 
r Xj1 I 



1/2 

1 13 3 

-^i(4)^>E^°(pq) *«u|(m) 

^odd 
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(17) 



For the sake of a simpler notation we used m = 1/2. As shown in Appendix B, this can be cast into the form 



4 pm=l/2i 
^1/2 ' 



= -EQ(p-q)*i;n|(w) 



1 

47T 



<t(12)--(<x(3)x«x(12)) 



p x q 

pq 



x m=1/2 >. 



(18) 



Of course, this relation is also valid for m = —1/2 

Finally, we tui 
coefficients yields 



Finally, we turn to the last part of Eq. (JjJJ, the expression for 4 D™= 1/2 . A gain inserting the Clebsch-Gordon 



4 n m=l/2\ 
^1/2 I 



1 13 3 

^l(4)~)E^(pq) *u 2 i|(m) 
~3 



/A 

id-)---; 

10 1 v 2 ; 2 2' 



(19) 



Due to the overall positive parity and the antisymmetry of the state with respect to the two neutrons the sum over I 
and A splits as 



ZA21* 



E ^" *«aij + E * H+2 21| + E ^+2 « */+2 I 21 1 

^odd ^odd ^odd 



(20) 
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Now, however, the different types of coupled spherical harmonics are more difficult to split into simple second order 
expressions and scalar functions. This is elaborated in the Appendix C with the result 

3#*(M) = y^mA l +B l (Y 2ll (p)+Y 2 ^q)) 
^ 2 r+ 2 (Pq) = y^(m)Ci + D l Y 2ll (q)+E l Y 2ll (p) 
^ai(Pq) = ^ 2 f(pq)Q + A^(p) + ^^(q). 



(21) 



The quantities A\ , . . . , E\ are analytically known scalar functions depending on the momenta p and q. They can be 
inferred from Appendix C and the first relevant ones are given below in Eas. l|27|) and (|28p. Using the decompositions 
of Eq. (J2U in Eq. J2I|) leads to 



E yfx *iAaii = yli (pq) E x < + r ^p) E ^ + *v(q) E ^> 



(22) 



/A 



with 



A; = A Z * H21 3 + C; (*n+221§ + *i+2Z2lf 



(23) 
(24) 
(25) 



For the representation of the spin states | (lh)%Ms) of Eq. I|19fl in terms of the spin state | x m ) we use an equivalent 
but modified form as shown above. The details are shown in Appendix D. After some algebra we arrive at 



4 n m=l/2, 
^1/2 / 



I /A 

2V 2tt 



<r(12)- po -(3)-p 1 
- -<x(12)-<t(3) 



<r(12) • q cr(3) • q 1 



er(12) • <r(3) 



: m > E^ 

lodd 

x m ) E^ 



13 11 

2 4tt V5M 



<t(12) • q er(3) • p + <r(12) • p <r(3) • q 
-|p. q< r(12).ir(3)l | X m > E^' 



(26) 



It is now the time to compare our results to the scalar expressions given in [l5j . We see that the first eight terms 
of Eqs. (2)-(7) in [3 are identical to the ones derived here (see Eqs. © , , UHl, and JUJl). While in the 

scalar functions multiplying the scalar operators are unknown, here they are explicitly provided in terms of the partial 
wave function c omp onents calculated e.g. in a Faddeev approach. We want to point out that the last expression in 
Eq. (7) of Ref. is redundant. By itself it is not antisymmetric under the exchange of particles 1 and 2 (the two 
neutrons). It has to be multiplied by a scalar function which is formed from odd orbital angular momenta I. Doing 
this, one arrives after some algebra at the following result: The L = 2 piece is already contained in the previous three 
terms of | 4 D™ 2 ). The L = 1 piece is identically zero, and the L = part cancels among the two terms given in the 
last expression in Eq. (7) of 15] . 



B. Normalization 



In the previous subsection we started from a partial wave decomposition of the 3N wave function ^ m (p, q) (Eq. 10). 
By the very construction the individual terms are manifestly orthogonal. In their new forms, as derived in the previous 
subsection, this is no longer obvious. However, since the new forms are identical reformulations of the original terms, 
it is simplest to go back to Eq. @ to verify orthogonality and normalization. Then one sees immediately that Si, 

2 Pi, and 4 -Pi given in Eqs. i|13|l . (| 1 f>|> and (|18J) are orthogonal to each other. The normalization for those three 
pieces is given as 

( 2 S 1/2 + 2 F 1/2 + 4 Pi/ 2 I 2 S 1/2 + 2 P 1/2 + 4 P 1/2 ) = E E [ d PP 2 [ dqq 2 yf XLsS (pq)- (27) 
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The last term in Eq. ©, | 4 -Di/2), is more intricate in the form of Eq. I|26|) . Of course one can also go back to Eq. 10. 
However, instead of doing this, we go back half way to Eq. (|19|l and insert the decomposition given in Eq. I|22ll . This 
leads to the three terms 



4 jym 



1/2/ 



Y 2 (P) | (l|)|> 



ym i (i|)|> 



lodd 



yhm i (i|)|> 



lo 

1 k' m 



(28) 



which are in unique correspondence to the three terms in Eq. I|26|) . The question of normalization and orthogonality 
requires knowledge of the analytically known coefficients Ai , . . . , E\ inside Vj , Wj and X[ . 
As shown in Appendix C, they are given as follows. If we keep only / — A = 1, then 



Ax = 1 
B x = 



(29) 



and the coefficients C, D, and E do not occur. If we allow in addition I — 3 and A = 3, but no higher values, then 

A 3 



~^(25P 2 (p-q) + ll) 



B 3 = 

Ci = 

Dx = - 

Ex = 



35 
127 



fi(p-q) 



A(p-q) 



(30) 



Note that the coefficients Cx, Dx, an d Ex occur together with 1 = 3. The evaluation of the terms I = 5 and higher 
can be found from the general formula given in Appendix C. 



If one keeps only the I = 1 partial wave function component, then | D" l ^ 2 ) reduces to the simple expression 



li) \i=\- 



^(pq) I 



*1121|(W), 



which is orthogonal to the previous states and normalized as 



' ^ n m I ^ r> m \ 

i LJ X/2 I U X/2I 



J _ i = J dpp 2 J dqq 2 *i 12 i|(p«)- 



(31) 



(32) 



If one includes the I = 3 or A = 3 partial wave function components, which are in fact tiny contributions, one arrives 



at 



1/2, 



1=1,3 



Y 2 (p) | (li)§) 



^ 2 (q) I (4)f> 



4tt 
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V5 



+ 



+ 



^i(pq) I (i~)~> 



*332lS(M) + *3121«(P?) P l(P • q) 



11 /7 



4tt 



*1121|(P9) + Y Y 2*3321|(m) 



(*132li(P9) + *312lf(P?)) + y V 2*83131 J (P9)-P2(P ' 4) 



(33) 
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Now the three terms are no longer orthogonal to each other, but of course an identical reformulation of 



1/2/ 
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y?i(pq) I 



y 3 2 3 (pq) I (i|)|> 
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2 j 2 7 
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as given in Eq. Q. Consequently, the state | 4 Dy 2 ) is normalized as 



4 nm I 4 7~>m 



( 4 D 



1/2 



1/2/ 



=1,3 



+*?aai*(P9) + *Li8(w) 



(34) 



(35) 



The direct verification of this result in form of Eq. (|26|l is straightforward but tedious. 

Summarizing this section, the sum of the expressions in Eqs. (jHJ, (|13(l . Ijl6|l . (|18|l . and (|26|l is the operator form of 
the 3N bound state in momentum space we were looking for. It has the form 
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* m (p,q) = X>(p,q)lxi>, 

i=l 



(36) 



where the \xi) are composed out of 8 scalar operators acting on the special 3N spin state |x m ), introduced in Eq. 
For the convenience of the reader, we list the states \xi) again. 



1X2) 

1X3) 
1X4) 
1X5) 
1X6) 
1X7) 
1X8) 



— <x(12) • er(3) I X m ) 
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2 i 



1 1 



pq 
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-j=- (i<r(12) p x q ( CT (3) x «r(12)) • (p X q)) | X r 

P x q 

pq 



-<r(12) • <r(3) 



-KI 


«r(12) - |K3) x 




V(12). p< r(3)-p 










r<r(12).q<r(3)-q 






q 2 


3 1 
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<r(12) qtr(3) 


2 v/5jp« 



-<r(12) • er(3) 



(37) 



Each term is composed of scalar operators consisting of spin operators and momentum vectors, applied on the pure 
spin state | x m )i which carries the overall quantum number J = 1/2. Furthermore, each of those terms in Eq. 13GI) 
includes scalar functions formed out of the two Jacobi momenta p and q. Their dependence on the standard partial 
wave function components has been explicitly worked out and will be investigated in the following Section III. Of 
course, exactly the same forms are valid in configuration space. In this case the Jacobi momenta p and q would have 
to be replaced by the corresponding conjugate configuration space Jacobi vectors. 
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III. THE SCALAR FUNCTIONS 



The operator form of the 3N bound state as given in Eq. I|36[) contains the scalar functions 4>i. They will be 
investigated now and are given according to Eqs. ©, i|13|) . (|18|l . and (|26[1 by 



*>i(p,q 
^(p,q 
^(p,q 
^(p,q 
^(p,q 
^(p,q 
Mp,q 
^(p,q 



i- 53 V^TTp^pq) *uooi(P9) 
i- 53 \/27+T p ( Pq ) * H01 1 (p?) 

^odd 

^ 5Z c Kp-q)*Hioi(w) 

^odd 
^odd 

^= 53 [Bi * lm s{pq) + Ei l+2 21f (pg) + D, * ;+2 , 21 g (pg) 

'odd 

— 53 (S z * H21f (pg) + Di ;+2 21f (pg) + E l *, +2 , 21| (pg) 

'odd 

4^ 5Z ( A < *«mf (w) + C| (*, (+2 2lf (pg) + * ;+2 , 21f (pg) 



4tT 



(38) 



These functions (f>i(p, q) depend on three variables, p, g, and p • q. They are determined by the partial wave 
components of the 3N bound state and analytically known coefficient functions. In Figs. [I] through |H we display 
4>\ (p, q, cos 9), 08(p,9,cos#), <^7(p, q, cos0),and 2 (p,g,cos#) for a fixed angle = 0. We see that the numerically 
largest function is </>i(p, g, cos 9), the other three ones shown are at least an order of magnitude smaller. While 
4>i (p, g, cos 9) has a simple, bell- like shape with the maximum at p = and g = 0, for 0g(p, g,cos0) the maximum is 
shifted to p w 0.2 and g ss 0.4 fm _1 . The reason for this is that (f>s(p, q, cos 9) does not contain s-wave contributions 
but instead includes tensor force couplings. The function 07 (p, g, cos^) is similar in shape to <ps(p, q, cos 9). Finally, 
cj>2 (p> <7> cos 9) also has its minimum shifted away from the origin. 

The dependence on the angle p • q is generally rather weak. In order to show the angular dependence explicitly, 
the function 0i (p, g, cos 9) is displayed in Fig. [S] for some fixed values of p and g as function of cos#. Similarly we 
show in Fig. the angular dependence of 4>%(j), g, cos#). The angular dependence of 07 (p, q,cos9) and 2 (p, g,cos#) 
is dominantly given by P± (cos 9), and thus not displayed. 

For the calculations presented in the following, we used a wave function based on the NN force AV18 0] in 
conjunction with the Urbana-IX three-nucleon force We show results for the 3 He. The 4>i functions for other 
interactions are qualitatively similar, especially, their relative importance is not changed. Tabulated functions for 
several force combinations are provided by the authors [2(J ■ 

In order to quantify the relative importance of the eight functions </>, (p, q) we consider the normalization of the 3N 
bound state 



tip 



I*' 



= 53 / d P rf q {Xi\Xi)$ + 2Pe 53 / dp dq ( X 



\Xd) 



(39) 



using the representation given in Eq. i|36|) . The numerical evaluation is straightforward, and the contributions of the 
different products of the </>i(p, q) (denoted as Ni 3 •) to the norm are listed in Table I. Clearly, the major contribution to 
the norm (91.42 %) is given by 4>\(p, q). The second largest contribution is already more than one order of magnitude 
smaller and is given by </>|(p,q). All other contributions are even smaller. 
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IV. APPLICATION 



As example for the application of the above derived operator form of the 3N bound state we consider the spin 
dependent momentum distribution of a neutron inside a polarized 3 He nucleus. This quantity is defined as 



iV(q) = <* m =5| ( 5( q _ qop )I(i + (To (3))|^ 



(40) 



It should be noted that in the 3N cm. system the Jacobi momentum q is the momentum of one nucleon, here nucleon 
3. Regarding the eight operator structures displayed in Eq. (|37|l one recognizes that only the following five different 
terms occur 

O x = 1 

2 = er(12)-er(3) 

3 = <r(3)-A=0 3 (A) 

4 = er(12)-B = 4 (B) 

5 = o-(12)-Co-(3)-D = 6 (C,D). (41) 

Here the vectors A, B, C, and D represent different momentum vectors. It is a straightforward exercise to evaluate 
once and for all the matrix elements 



The nonvanishing ones are listed below: 



(42) 



(Old) 


= 1 




(Oi0 3 ) 


- A Q 




(o 2 o 2 ) 


= 3 




(0204) 


= B Q 




(o 2 o 5 ) 


= C 


D + i(C x D) 


(0301) 


- A 1 

— ^0 




(o' 3 o 3 ) 


= A' 


A + i(A' x A) 


(0402) 


= B> 




(04O4) 


= B' 


B 


(Q405) 


= (B' 


■C)D 


(o' 5 o 2 ) 


- C 


D'-i(C' x D')o 


(0504) 


= (B 


C')^o 


(o' 5 o 5 ) 


= (C 


• C)[D' • D + z(D' x D) 



(43) 

For the evaluation of the specific expectation value N(q) considered in Eq. (|40|l we need in addition matrix elements 
of the form 



(0>o(3)0 J )Mx m |0^ (3)0,|x m ). 



(44) 



The resulting, nonvanishing matrix elements are listed in Appendix E. With this, the specific operator O for evaluating 
N(q) from Eq. (|40J) has spin matrix elements given as 



(0[-(l + a Q (3))0 3 ), 



(45) 



and the nonvanishing matrix elements arc also listed in Appendix E. Next one expresses the states \x%) in terms of 
the 5 operators Oi, 



\xi) = Y, A M^)\x m ), 



(46) 
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where denotes the arguments of the operators Oj which varies with \Xi)- One explicitly obtains 



lx 2 ) - ^o 2 1 x m ) 



1X3) 



3 1 Q 3 (pxq) 
2 i pq 



x m ) 



1 1 



|X4) = -=-(i0 4 (pxq)-0 5 (q,p) + 5 (p,q)) | X m ) 

V2 w 



1X5) 
1X6) 

\Xr) 
Ixs) 



1 ^ 4 (p x q) _ i Q 5 (q,p) - Q 5 (p,q) 

2 ------ 



p<7 

Q 5 (P,P) i n 
Q 5 (q,q) l n ' 



pq 

x m ) 
x m ) 



x m ) 



3 1 1 

2 VEm 



5 (q, P) + O s (p, q) - -p • q0 2 



x m )- 



(47) 



Using the expectation values listed in Appendix E, one can determine the matrix elements (xi\0\xj) m a straightfor- 
ward fashion. As example we give 



(X4|0|xe) = ^r — (^(P x q) - ^P • q + p qo)- 
2 pq 3 p z 



(48) 



Here the index denotes a spherical component. The above expression nicely exhibits the analytic angular dependence. 
The final step in obtaining the momentum distribution A(q) is then to write 



N(d) = J dp (* m =*(pq)|6|# m =*(pq)) 

dp (5>?<X*|0|Xi) +2^^ J Re( Xi |6| Xj ) 



(49) 



An inspection of all analytically given spin matrix elements reveals that only four types of angular integrations 
occur. These are 



/ 



dp f(p ■ q) < 



1 

Pq 

vl 



(50) 



I (P x q)o 



Because of the rotational invariance around the quantization axis (z-axis) one can choose the vector q to be in the 
x-z plane. Then it is most convenient to rotate the z-axis into the direction of q by the angle 9 = 9 q . The result is 



1 



1 



/ dp/(p-q) H = ^ ftef{x) \ P p 



(P x q) 2 , 



cos Ox 

2 (cos 2 6x 2 + ± sin 2 (9(1 -x 2 )) 



(51) 



\p\l-x 2 ) 



Collecting all the results one ends with the final expression for A^(q) = N(q, 9) given by 



dx 



N(q,6) = 2tt / dpp 2 
Jo 

xjy?(p, q, x) + ^4>l{p,q,x) 



+4> 2 (p,q,x)<p 6 (p,q,x)-^=(3cos 2 9 - l)(ix 2 - 1) 
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\[2 

+(f> 2 (p,q,x)(t) 7 (p,q,x) — (3cos 2 6>- 1) 
2 

+02 (p, q,x)(j)s(p, =(3cos 2 6» - l)x 

vlo 



+ ^(p ;(? , x )^(l_ C o S 2 0)(l-:r 2 ) 



+0 2 (p, g ,x)-(3 + cos 2 #)(l-z 2 ) 



+(f) 4 (p,q,x)(t) 5 (p,q,x)-^=(coa 2 9 - 3)(1 - x 2 ) 
V3 

+04 (p, 9, x)<f> 6 (p, q, x ) — ( 3 cos2 - 1)(1 - ar 2 )^ 
3 

+^Aip,q,x)^{jp,q,x)— = (3cos 2 6> - 1)(1 - a; 2 ) 
zy ID 



+0 2 (p, g ,x)i(cos 2 + 9)(l- 2 ; 2 ) 

+0 5 (p,9,x)0 6 (p, ( j,x)iy|'(3cos 2 6'- 1)(1 -x 2 )a; 
3 

+05 (p 5 9, z)0 8 (p, 9, a;)— =(3 cos 2 6* — 1)(1 — a; 2 ) 
4V5 

+0§(p, 9, as) ^[4 + (3 cos 2 9 - l)(ix 2 - 1)] 

+06 (p, 9, :r)0 7 (p,g,a:)- (3 cos 2 9 + l)(3x 2 - 1) 
6 

+06 (p, 9, z)0 8 (p, 9, s)t4=[8 + (3 cos 2 - l)(3a; 2 + l)]x 
2 V 30 



+0 2 -(p,9,z)-(3cos 2 + l) 
6 

j~~2~ 

+<fa(p,q,x)<ps(p,q,x)\ — (3cos 2 6>+ l)x 
V 15 

+0 2 (p, q, i) -L [48 + (15 cos 2 9-1) {bx 2 + 3)] } . (52) 

Though the angular dependence for the direction of the nucleon momentum q in relation to the quantization axis is 
analytically given, the full expression is quite lengthy. However the dependence on the angle 9 is quite simple, namely 
N(q, 9) = a(q) + b(q) cos 2 9. Moreover most of the contributions are numerically insignificant, as illustrated in the 
previous section. Therefore, we only keep 0i, 0s, 07, 02- It turns out that for the specific quantity N(q,9) only these 
components visibly contribute. Therefore, in this case the lengthy expression of Eq. (|52[1 shrinks to the few leading 
terms 



/>oo />1 

AT (9,0) = 2tt y dpp 2 J dx 

x{0i(p,9,a;) + ^02 (p, 9- x ). 



+<f) 2 (p,q,x)<pr(p, q, x) — (3 cos 2 9 - 1) 
2 

+02 (p, g,x)0 8 (p, 9, x)—=(3cos 2 9 - l)x 
V 15 



+0 2 (p,9,x)i(3cos 2 + l) 
/ 2 

+07 (p, g,x)0 8 (p, 9, a;)v — (3cos 2 6» + l)x 
V 15 



+08 (P. 9, x) 7^ [48 + (15 cos 2 - l)(5x 2 + 3)]} (53) 
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The numerical results for N(q) are displayed as function of q in Fig. for the fixed angle 8 = (i.e. the nucleon 
momentum is parallel to the quantization axis) and in Fig.|^|for 8 = 90°. We compare the full result given in Eq. ijol'l) 
to various truncated sums. The solid line in Figs. and corresponds to the full calculation using Eq. I)52[l . The 
most simple approximation would be to consider only the first term in Eq. Q52f) or Eq. (|53[l , namely <j>\ , given by the 
dashed line. We see that this simple term alone already gives a good representation of N(q, 8) up to about q s» 1 fm _1 . 
The dip around 2 fm _1 is mostly filled in by adding the term containing 0|, shown as dash-dotted line. When adding 
terms containing <f>>j, represented by the dotted line, one is very close to the full result. Adding the terms containing 
cj>2, i-e. calculating the expression given in Eq. (|53[1 . shows that all other terms in Eq. I|52|l are insignificant. 

We expect that also for other observables the operator form of the 3N bound state will be useful. It should provide 
an easy access to the wave function without the need of having access to a modern triton code. 



V. SUMMARY 

An old idea by Gerjuoy and Schwinger [l^ has been revived to present the 3N bound state in operator form. 
This form analytically exhibits the dominant angular and spin dependence of the wave function in form of scalar 
operators formed out of momentum and spin vectors, which are applied on a pure spin 1/2 state. Each such operator 
is accompanied by a scalar function depending on the magnitudes of the two Jacobi momenta and the angle between 
them. We established the connection of this form with the standard partial wave decomposition. This connection 
provided the explicit form of the scalar functions in terms of partial wave function components. The key point in the 
derivation was, to extract from an infinite sum of partial wave expressions the operator form and the accompanying 
scalar functions. The presented operator form of the 3N wave function is independent of the applied NN and 3N force. 

We illustrated the application of this new form of the 3N bound state wave function by calculating the spin 
dependent single nucleon momentum distribution in a polarized 3N bound state. It turned out that for this quantity 
only 4 parts out of the total number of 8 parts forming the 3N bound state were needed to achieve a sufficiently 
accurate representation. Several sets of spin matrix elements depending on the Jacobi momentum vectors, which have 
to be calculated only once, have been evaluated. They will also be needed in other applications. 

We expect that this operator form allows an easy access to the 3N bound state. The 8 scalar functions carr yin g 
the specific dynamical information have been tabulated on a sufficiently fine grid and can be downloaded from [20j. 
There are sets of scalar functions for various modern NN and 3N force combinations. 
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APPENDIX A: RELATION FOR THE COUPLED SPHERICAL HARMONICS 

Here the relation of Eq. (|llfl for the coupled spherical harmonics 3^/'(pq) will be verified. We consider 
3- «° (pq)y 'n (pq) i which can be rewritten by standard techniques (see e.g. ^(|) as 

^ (pq)<(M) = ^-V2lTT^(-r +i c(/ia,oo)c(m,oo){ \\\ \y%m- (ai) 

71 ab ^ J 

The sum over a and 6 will give four terms. After inserting the explicit expression for the Clebsch-Gordon coefficients 
and the 6-j symbol, one arrives at 



■uOO/~v»;i/v-M 3 1 1 / (I + !)(/ + 2) ^ 1(1 + 1) ^ \ , 40 , 



Using the relation 



^°(pq) = (-) Z ^^fl(p-q) (A3) 



the previous Eq. (| A2|) can be rewritten as 



2Z + 1 

(pq) = V^TT) 



^ i i-)'+V|(2'-i)^-i(p-q)yif(pq) + ^/ (i 2 ) )a 3 V -V,<pqi | ■ (A4) 



This is a recursive formula and leads to the relation given in Eq. (II ip - In practice only low orbital angular momenta 
occur and one easily works out the lowest terms as 

ci(pq) = 1 

C2(pq) = -\/5Pi(p-q) 

c 3 (pq) = ^y|(5P 2 ( P -q) + l). (A5) 



APPENDIX B: VERIFICATION OF THE RELATION GIVEN IN EQ. (tl8l) 

The first step is to generate the spin states of Eq. I|17|) from the states | x m ) given in Eq. JSJ. It is easy to see that 

\(l^) = a(12) + \ X m= h, (Bl) 

where cr+(12) is the spherical component of the spin operator given in Eq. ©. 
This form can be rewritten as 



= |(o- + (12) - |(5(3) x a(12)) ) | X ro=i ) (B2) 



using the spherical component of the vector product. 

Next, we express the second spin state in Eq. I|17|) with the help of the lowering operator as 

l(l^)~) = (S x ~iS y ) I 



= -L( < 7_(l)+ < 7_(2)+<7_(3))i(<7 + (l)-<7 + (2)) | X m =^> 



Finally, starting from 



i=i(2(a (l)-a (2) + a_(3)( ( T + (l)-a + (2))) \ X m =?) 
-^1(2(0-0(1) - ffo(2)) S{2)) x c?(3))J | X m ^> 

A( CTo (i2)_ 1(^(3) x^(12)) ) |x m=i >- 



inserting the form (B.3) and reshuffling leads to 



1 \ } = ~ X ° ;(12)) - ] 1 X 



ra — - r 



Now we use the relation (2.f0) and the property 



•vW-m • 3 1 (P x <?V 



With this one arrives directly at Eq. (jf 8|) . 



APPENDIX C: VERIFICATION OF THE RELATION GIVEN IN EQ. pljl 

With standard recoupling techniques one finds 



v oo- Aw2mr6fil _ 3 f / / q+l)(/ + 2)(2Z + 5) 2m 
% (pq)y n (pq) - W 30(2/ + 3) 



—5— + ^-h+j + y (2z-i). 30 ^- 1 ><- 1 J ' 

which leads to the recursive relation 



,,2m _ / 30(2/ + 1) J 4tt , oo 2m 

% " t/j(l + l)(2I + 3) y^-^wi-^n 



5 y (2Z - 3)30 ^- 2J - 2 

Similarly, starting from 



?™(f>h)Y fnl - 1 L/ 3 ([+1)0 + 2) 2m 

(pq) r ^(p) - 7^1 V2 V(2TTT)(2TT3) fe 



' Kj + 1) V 2m 1 /3 / *(j ~ 1) V 2m 

(21 - l)(2l + 3) " V 2 V (21 - l)(2l + 1) l ~ 21 



one finds the additional recursive relations 



V^rAtfl / 2 / (2^ + l)(2Z + 3) /T- v o 0y _ 



' *(j + 1) y2m /3 / Kj - 1) V 2m 

(2Z - 1)(2J + 3) " V 2 V (2Z - l)(2l + 1) ^ l - 2 ' 1 
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and similarly 



3^2 (pq) 



f (2/ + l)(2l + 3) 
(I + !)(! + 2) 



47r^°y 2m (q) 



(2Z — 1)(2Z + 3) 



yu -\/o 



1(1 -1) 



(2Z- 1)(2/ + 1) 



iJ-2 



(C5) 



Inserting these equations into each other yields the relations given in Eqs. (51) ■ For the calculation of a 3N bound 
state, the scalar functions Ai to E\ are in practice only needed for small values of I. In our context we only need odd 
values of I. The two lowest cases are given here. 

If only I = 1 is kept, then one trivially has A\ = 1 and B\ = 0, and all other terms are absent. If only I = 1 and 
I = 3 are kept, then one obtains from Eqs. (C.4) and (C.5) 



Jzm / ~ ~ \ / -\?2r, 

31 (Pq) = Vg^l! 



2 m 



yiT(pq) 



2m 



_5 

47T 

~5 

47T 



PiY 2m (p) 
PiY 2 m(q). 



Furthermore, Eq. (C.2) yields 



34 m (pq) 



'^{yfnm^Pi 



6 (yfr+yir) 



and after insertion of the results (C.6) we get 



%3 (PQ) = 2 ^( 25p 2(pq) + 11) - 

Thus, one obtains in the end the coefficients given in Eq. I|30(l . 



35 

12^' 



Pi(pq) U2m(p) + F 2m (q) 



(C6) 



(C7) 



(C8) 



APPENDIX D: THE OPERATORS OF EQ. lt26l) 

The starting points for the derivation are Eqs. (|19fl and l|22|l. According to Eqs. I|B2(1 . i|B3|) . and (|B5|> these lead to 
I 4 ^ =i > = -^y^(a(12)-i(a(3)xa(12))) + | X m ^> 



]T (^V^pq)^ + r 2 -i(p)^ + iww 



lodd 

T 2 



X] (^n (pOT + ^20 (p)Vj + y 2 o(q)W( 



lodd 



-—(a(12)- l -(a(3)x*(12))_\ X m = L *) 
(^(pq)X, + r 2 i(p)Fi + y 2 i(q)W,) 

y|<7_(3)<7_(i2) | x m =5) £ (^(pq)^ + r 22 (p)v^ + y 22 (q)Wi). (Dl) 



lodd 

[2 / I 
V 5 



lodd 



As an example let us regard the terms in Y^p), which have the well known representation in terms of spherical 
components of p: 

v r*\ 1 I™ 
y 2 _i(p = — pap- 

2p2 V 7T 
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Y 2 i(p) 



2p 2 
1 / 1 r, 



2p^ V 7T 

For the convenience of the reader we also provide the relations 

1 3 1 



yn°(pq) 
^n(pq) 



V2 47T I p || q 
3 



(PoQ- +P-qo) 



111, 

-7={p+i- +p-q+) 



4tt | p || q | \y/Q 
1 3 1 



V2 4tt | p || q 

3 p+q+ 



{p + q +p q+) 



(D2) 



(D3) 



4tt | p || q | 

If one now inserts (D.2) into (D.l) and looks only into the terms with Vi, one can easily combine the expressions to 



4DT 



v L 2 V 2tt ^ 

lodd 



<r(12) • p g(3) • p 1 

^2 3^(12) -a(3) 



(D4) 



The term in X\ is somewhat more tedious. 



APPENDIX E: THE NONVANISHING MATRIXELEMENTS OF EQS. (gij| AND (1551) 



The nonvanishing matrix elements of Eq. 14411 are given by 



(01(70(3)0!) = 1 

(Ok7o(3)0 3 ) = A) 

(O 2( 7o(3)0 2 ) = -1 

(O 2 (7o(3)04> = B 

(O 2( t (3)O 5 > = 2O A) - C • D - i(C x D) 

(O^ (3)Oi) = A' 

(O' 3 a (3)O 3 ) = 2A' A - A' ■ A - i(A' x A) 

(O 4 a (3)O 2 ) = B' Q 

(O 4 a (3)O 4 ) = B'-B 

(O 4 a (3)O 5 ) = (B'-C)A, 

(O' 5 a (3)O 2 ) = 2C' Q D' Q C D' + i(C x D')o 

(O^ (3)O 4 ) = (B-C')D' 

(O' 5 a (3)O 5 ) = (C'-C)[2jy o i) -D'-D-»(D'xD)o]. (El) 
The nonvanishing matrix elements of Eq. Q45J) are given by 



(OiOOi) = 1 

(OiOOg) = A 

(o 2 6o 2 > = i 

(0 2 00 4 > = B 
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\ ■■■> i 


= C D 


(O'.OOi) 


U 


(o;6o 3 ) 


U u 


(o' 4 do 2 ) 


= 

U 


{o' A dOi) 


= B' B 


(o'.dOc,) 


= fB' • C)Ai 


(0^00 2 ) 


= <M 




- (B-C'K 




- (C'-C)D' D 



(E2) 
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TABLE I: The contributions Nij from 4>i(p, <l)(j>j(p, q) to the total normalization of the 3N state. 



ij 


11 


22 


33 


44 


55 


66 


67 


68 


77 


78 


88 


Nij [%] 


91.42 


0.76 


0.02 


0.02 


0.02 


0.13 


0.11 


0.25 


0.98 


2.44 


4.35 
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FIG. 1: The dependence of the function <j>\(p,q,x = 1) on p and q. For better visibility also contour lines are given. 



0.06 
0.04 




FIG. 2: The dependence of the function 4>s(p, q, x = 1) on p and q. For better visibility also contour lines are given. 





FIG. 4: The dependence of the function <j>2(p,q,x = 1) on p and q. For better visibility also contour lines are given. 
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(x 50)p = 1.15fm" 1 ,q 



: 1.14 fm" 
: 0.25 fm" 
: 1.14 fm 



COS 
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FIG. 5: The angle dependence of (j>i(p,q,cos8) for different pairs of momenta p and p as indicated in the figure. Note the 
multiplicative factors in three cases. 
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FIG. 6: The angle dependence of <f>s(p,q, cos 8) for different pairs of momenta p and p as indicated in the figure, 
multiplicative factors in two cases. 
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1e+1 




1e-6 1 1 1 1 1 1 

1 2 3 4 5 



q [ f m ] 

FIG. 7: The momentum distribution N(q, 9 = 0). The solid line represents the calculation containing all terms in Eq. 1521 , 
The dashed line displays the result based on the first term only, the dash-dotted the one containing in addition the contribution 
of </>&. For the dotted line contributions containing <j>f are added. The long dashed-line contains in addition contributions 
including <j>2 and corresponds to Eq. H53H . 
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FIG. 8: Same as Fig. 7, but for N(q,0 = 90). 



